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Abstract 

I discuss three subjects in thermal field theory: why in SU (N) gauge 
theories the Z(N) symmetry is broken at high (instead of low) tempera- 
ture, the possible singularity structure of gauge variant propagators, and 
the problem of how to compute the viscosity from the Kubo formula. 

Based on an invited talk presented at the Workshop on Hot QCD held at the 
Winnipeg Institute for Theoretical Physics, Winnipeg, July 1992. 



I take this opportunity to make three comments about field theories at nonzero 
temperature, especially as applies to hot QCD. The first is trivial; the second, base 
speculation, and the third, a technical problem. 

1. Why the Z(N) symmetry is broken at high temperature 

Consider an SU(N) gauge theory without dynamical fermions. At nonzero tem- 
perature the order parameter for this system is the Polyakov line [ffl], 



L(x) = tr exp ^ig J A (x)dr 



V refers to path ordering, g is the coupling constant, and (3 = 1/T, with T the 
temperature. The Polyakov line is invariant under gauge transformations which are 
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strictly periodic in the euclidean time (3. There are also gauge transformations which 
are periodic only up to an element of the center of the group 0; the aperiodicity is 
a global Z(N) phase, = exp(i<pj), 4>j = 2nj/N with j = 0...(N — 1). The Polyakov 
line transforms under these aperiodic gauge transformations as L(x) — > exp(i<pj)L(x). 
The Polyakov line has a nontrivial vacuum expectation value at high temperature, 
choosing one of N equivalent vacua, and vanishes below the transition temperature. 
This behavior is exactly opposite to that of a typical spin system, where the symmetry 
is broken at low temperature, and restored at high temperature. Here I give a simple 
argument why. 

In a spin system with discrete symmetry, the simplest argument for spontaneous 
symmetry breaking at low temperature is due to Peierls [|3|. Consider an Ising model, 
with a spin = ± on each site, and a coupling J between adjacent sites, so the hamilto- 
nian = — J °i ' Start with a perfectly ordered state, where all spins = +. Now 
insert a domain of width ~ R into this state, where all spins = — inside the domain. 
Because the hamiltonian vanishes when the spins are aligned, the only energy is from 
the interface of this domain, = J A, with A ~ R d ~ x the area (d is the number of 
spatial dimensions). This domain contributes to the partition function as 

Z ~ V e- JA ' T . (1) 

The factor of volume in the measure, V ~ R d , is from the entropy of moving the 
domain around. Peierl's argument is that in more than one dimension, at low tem- 
peratures the factor of R d_1 in the exponential always wins over the R d in the measure: 
domains are exponentially suppressed, and the symmetry is spontaneously broken. In 
this view, it is not apparent how the symmetry is restored at a finite temperature; 
this happens because when domains are common they interact with one another, and 
this lowers their total energy [[| . 

Z(N) interfaces can be constructed in hot gauge theories. In three spatial dimen- 
sions the interface tension between two Z(N) domains is @] 

T 4(N - 1)tt 2 T 3 
3V3N g 

The detailed origin of the interface tension is a story in itself. This result is com- 
puted semiclassically, and so is valid at high temperature, where the running coupling 
is small. The interface tension is the stationary point of an effective action, and rep- 
resents a balance between classical and (certain) one loop effects. Because of this 
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balance, J is proportional not to 1/g 2 , as is standard for instanton type solutions, 
but to 1/g. 

For my purposes here, however, the powers of g are inconsequential - - what 
matters is that the mass dimension is set entirely by powers of temperature, and 
not (say) by the renormalization mass scale. It is then obvious why the symmetry 
is broken at high instead of low temperature: the effective spin spin coupling is not 
constant (as for the Ising model), but increases with temperature, growing as ~ T 3 at 
large T. The contribution of a domain to the partition function is Z ~ exp(—#T 2 R 2 ) 
at large R and T: interfaces are rare at high temperature, and become common only 
as the temperature decreases. 

What happens at large N, when g 2 N is held fixed? The usual free energy of a 
gas of gluons is proportional to iV 2 , so it is natural to expect the same for the spin 
spin coupling J. Instead, one sees from (|2]) that if g 2 N is held fixed at large N, that 
J ~ N and not N 2 . But remember that the Z(N) vacua are for phases <pj = 2irj/N, 
so if two Z(N) domains differ by a finite amount in j, the change in phase is only 
of order 1/N. For two Z(N) domains to differ by a finite phase, the difference in j 
must be of order N. The total energy for such an interface is of order NJ, which is 
of order N 2 , as expected. 

There is an amusing corollary to this picture. Consider the above argument at low 
temperatures. If the interface tension is nonzero as T — > 0, then Peierl's argument 
implies symmetry breaking about zero temperature. Assuming that the symmetry 
is unbroken about zero temperature, then for the picture of Z(N) domains to apply 
it is necessary for the interface tension to vanish at least linearly with temperature, 
J ~ T as T -> 0. 

I end this section with a matter of definition. I propose that the vacuum ex- 
pectation value of the Wilson line is related to the free energy of a test quark, F q , 
as 

(L(x)) = e - F « /T+l ^ . (3) 

The phase factor depends upon which Z(N) domain the theory falls into. This 
relation is normally written without the phase factor . But then it is only consistent 
in the trivial phase, <pj = 0; otherwise, @ implies that when j ^ 0, the free energy 
is complex, which is absurd. 

The phase factor in (||) can be understood from the definition of the underlying 
partition function. To summarize: the Polyakov line represents not just the "mass" 
of the test particle (which at finite temperature is the free energy), but more gen- 
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erally, the propagator, L ~ $(/?)V>(0). Now the gauge field, which resides in the 
adjoint representation, is invariant under the aperiodic gauge transformations men- 
tioned previously. Fields in the fundamental representation, however, do change: while 
ip(0) — ► ^(0), ${0) — > exp(i(j)j)ip(P). This gives exactly the transformation expected 
for the Wilson line. More generally, one can show that the transformations of the 
Polyakov line under Z(N) transformations reflect the invariance of the states of the 
(purely gluonic) theory under global Z(N) gauge rotations 0. 

2. Singularity structure of gauge variant propagators 

In this section I conjecture what the possible singularity structure of gauge variant 
propagators are like in covariant gauges at nonzero temperature. The motivation is 
recent work by Baier, Kunstatter, and Schiff |J, who showed that special care must 
be taken in covariant gauges. 

I first give a simple argument for the singularity structure of gauge variant prop- 
agators at zero temperature. Consider QED in covariant gauge, with gauge fixing 
parameter £. Ignoring Dirac structure, at one loop order the fermion self energy for 
a fermion of mass m has the form 

£ lioop (cu,A;)~ ce 2 (u-E k ) H ^ = (u - E k ) ce 2 In ( ) . (4) 

Ju-E k K \U — E k / 

The external energy is assumed to be near the mass shell, u ~ E k , and c = (3— £)/167r. 
The logarithmic divergence is usually cut off by introducing a photon mass or some 
such. Instead, let me leave it as is, and assume that the effects of all higher order is 
simply to sum up the lowest order term into an exponential. Then the logarithmic 
divergence at one loop order becomes a branch point: 

A~V, k) = uj - E k + Z lloop (u, k) + ... ~ (w - E k ) 1+ce2 . (5) 

While the above is just a guess, the result is correct: the correct singularity structure 
for the fermion propagator in QED is a branch point, where the strength of the 
singularity is gauge dependent 0. 

At nonzero temperature I simply assume that the form is similar to that at zero 
temperature, with the addition of the appropriate factors for a thermal distribution: 

Y} loop {uj,k)~{u-E k )tce 2 ^(l + 2n(A;)). (6) 

Jw~E k k 

Here n(k) = l/(exp(k/T) — 1) is the Bose-Einstein statistical distribution function, 
and c is a computable, nonzero constant. Note that at nonzero temperature the 
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gauge dependent term in (|6|) is strictly proportional to £, and vanishes if £ does. This 
is because at nonzero temperature, the physical longitudinal and transverse modes 
acquire thermal "masses" which themselves cut off the infrared divergences near the 
mass shell. The sole exception to this are gauge dependent terms in covariant gauges; 
there the poles, massless at tree level, remain massless to any order in perturbation 
theory. Also, I assume the upper limit on any infrared divergence is now set by the 
temperature. It is then easy to see that for nonzero temperature, since the Bose- 
Einstein distribution function behaves as n(k) ~ T/k at small k, that the integral for 
the self energy now has a linear divergence, 

E 1 {u-E h )£cf? (-^=r) ■ (7) 

\uj - E k J 

This integral is precisely that studied by Baier, Kunstatter, and Schiff || : the factors 
of us — E k appear to cancel, giving a gauge dependent contribution even on the mass 
shell. Baier, Kunstatter, and Schiff showed that this term has an imaginary as well 
as a real part, and so appears to give a gauge dependent contribution to the damping 
rate — which Braaten and I claimed was gauge invariant ||. 

Rebhan was the first to point out that the above integral is infrared singular [BJ; 
if computed with an infrared regulator, the term vanishes on the mass shell. This is 
simply because at one loop order there is a delicate balance between uj — E k upstairs 
against one over the same factor from the integral. Any infrared regulator will tend 
to smooth out the integral, so the uj — E k upstairs dominates, and Y} loop vanishes on 
the mass shell, uj = E k . 

My purpose here is merely to note that this is very special to one loop order and 
fails at higher order. Indeed, let us assume — as is true at zero temperature — that I 
can sum up the most infrared singular terms by exponentiation of the one loop result, 
to give 

A-\oj, k) ~ {uj - E k ) exp (lt~ce 2 T —^\ . (8) 

V uj — E k ) 

Because of the power like infrared singularities at nonzero temperature, instead of a 
branch point singularity there appears to be an essential singularity. 

I do not have a grand moral from this calculation. Probably for reasons of prej- 
udice as much as anything else, I believe that the position of the singularity in a 
gauge variant propagator is a gauge invariant quantity. Certainly the general proof 
of Kobes, Kunstatter, and Rebhan JlO| does not depend on the detailed nature of the 
singularity in the inverse propagator: only that it vanishes on the mass shell. But the 
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tricks that work at one loop order probably will not be good enough beyond that. A 
more general methodology is needed. 



3. The problem of viscosity 

It's really quite surprising that it's so difficult to compute damping rates in hot 
gauge theories. Inevitably, this has bred a feeling of ennui and despair: after all, what 
are the damping rates good for? How do they affect measurable quantities? 

One class of quantities which depend crucially on the damping rates are the trans- 
port coefficients | 11| , |12| . In this section I repeat calculations of Ilyin et al. (11 
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order to emphasize an apparent contradiction. Consider the Kubo formula for the 
shear viscosity, rj; the bulk viscosity, (, can be computed from rj. The Kubo formula 



for rj is [1 1 



1 1 

n = - lim — Im D(-iE + + ,0) . (9) 

D(p°,p) is the Fourier transform of a two point function of the stress energy tensor 
T ij : D(t, x) = < Tr ij (r, x) 7^(0, 0) >, where n ij = T ij -T kk 5 ij /3. The discontinuity 
is computed after analytic continuation from euclidean p° to p° = —iE + + . 

At lowest order the two point function of T u is given by the exchange of two 
transverse gluons, 

D{p°,0) = —TrA t {k ,k)A t {p -k°,k)(7k' 1 - 10k 2 k°{p°-k°) + 7(A; ) 2 (p -A; ) 2 

^ 3 (10) 
Here k = \k\, Tr = TJ2k° I d 3 k/(2ir) 3 is the integral over the loop momentum, and 
At is the propagator for transverse gluons. The four powers of momenta arise because 
each T u brings in two derivatives of the gauge field. Since the stress energy tensor is 
gauge variant, gauge variant modes do not contribute to flTUp. The contribution of the 
plasmon modes is neglected. This is allowed because momenta of order T dominate, 
and in this range the plasmon modes are negligible (their residue is exponentially 
small) . 

Using the spectal representation of the transverse propagator, 

rl/T n p+oo 

A t (k°,k)= dre lkT dcu p t (cu,k)(l+n(uj))e- UJT , (11) 

JO J-oo 

where p t (u,k) is the spectral density for transverse gluons, and n(uS) the Bose- 
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Einstein statistical distribution function, 



C roo r+oo / 



(12) 



The spectral density is assumed to be a Breit-Wigner form, 



r h) 1 f 7f(fc) 7f(fc) \ 

j "" 2ttA;v(w-A;) 2 + 7^) 2 + {u + kf + lt (kf ) ' 1 j 

where 7t(fc) is the damping rate for a transverse gluon of momentum k. Then 

C= - t^t r dk nik ^ + ^ = !r ? • < 14 > 

The contribution of quarks can be computed similarly. For three flavors of massless 
quarks, the total result in QCD is 

ltap = ^ /32 63\ 

225 7+7 ' 

In this expression 7 t and 7 + are the damping rates for gluons and quarks with mo- 
menta of order T, which in that range are independent of momenta. Computation 
shows that these damping rates are of order g 2 T, so in all 

HOOP _ ^ _ QgN 

9 

The exact value for r] lloop can be read off by using (|15f) and the damping rates in ref. 
JT3J . What I am interested in here are the powers of g. 

At first sight it might seem odd that the shear viscosity diverges in the limit of 
a free theory, g — > 0. In fact this is a familiar feature of the shear viscosity: it is 
infinite for an ideal gas because if an ideal gas doesn't start out precisely in thermal 
equilibrium, without interactions there is nothing to drive it there. 

The problem of viscosity is the following. The stress viscosity can also be calcu- 
lated using kinetic theory [I2|. As is standard in kinetic theory, the result for rj is one 
over the cross section. At lowest order the amplitude for the scattering of two gluons 
going into two gluons is of order g 2 , and so 

T 3 

^transport — , (17) 

Much work has gone into computing the subleading logarithms in (|17D ; Baym et al. 
12| pointed out that the scale is set not by magnetic mass (as first thought) but by 



7 



Debye frequencies, of order gT. But forget the logarithms — even the powers of g 
don't match up between (|1^) and (fPTI). 

This discrepancy is special to hot gauge theories. Consider a scalar field theory 
with four point coupling A. The damping rate first enters at two loop order, 7 ~ A 2 , 
so there one expects that both Kubo and transport theory give rj ~ 1 / A 2 . 

It was pointed out to me first by Gordon Baym, and later by Eric Braaten, that 
rj lloop is incomplete. This is also discussed recently by Jeon [fj|. The above calcula- 
tion represents the resummation of an infinite set of diagrams, which are self energy 
corrections on two hard lines in a loop. Besides this diagram, there are also ladder 
diagrams, in which soft gluons are emitted, all parallel to one another, from the top 
to the bottom hard line in the loop. One can easily show by power counting that 
these diagrams are singular as 7 — > 0, and so contribute to the stress viscosity. 

If the kinetic theory calculation is correction, then the ladder diagrams must 
dominate, and give 77 ~ l/f? 4 - Then r] lloop would be down by g 2 . Using the Kubo 
formula would just be a dumb way of doing things. 

Otherwise: the result is 77 ~ l/g 2 , with both ladder diagrams and rj lloo P contribut- 
ing at the same order; then there is more to the kinetic theory of ultrarelativistic 
plasmas than first meets the eye. 

I suspect that the former is most probably true. But in either case, one can NOT 
conclude that the damping rates are irrelevant; they are certainly part of the answer 
for the transport coefficients. 

Lebedev and Smilga |L5| have studied the contribution of the damping rate to the 
photon self energy in hot QED. They find that in the sum of the hard loop, similar to 
rj lloo P^ and ladder diagrams, that the damping rate cancels. While the damping rate 
cancels in the photon self energy, this does not imply that it does so in all quantities. 
Certainly there is a stress viscosity, computable (at least in principle) from the Kubo 
formula. 

I conclude this section by showing that once one has the stress viscosity in hand, 
the bulk viscosity follows easily in a surprising way. The standard relation between 
the bulk and stress viscosities is 

Here p is the pressure and e the energy of the gas at a temperature T; calculation at 
next to leading order shows 

p = T\ ai - a 2 g 2 + ...), (19) 
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where a\ and 02 are constants which depend upon the number of colors and flavors 
||16|| . The energy is given from the pressure by e = T 2 d(p/T)/dT. If the coupling 
constant is considered as a fixed parameter, then it is easy to see that even with the 
term in (|19|), that as for an ideal gas dp/de =1/3, and the bulk viscosity vanishes. 

Of course the coupling is not a fixed parameter: it "runs" , and thereby develops 
a nontrivial dependence on temperature. With c\ the lowest order coefficient of the 
/3-function, dg/dln{T) = —cig 3 (ci = (11 — 2iV//3)/(167r 2 ) for Nf massless flavors in 
SU(3) color), 



The number refers to three massless flavors in £77(3) color. 

The appearance of the lowest order coefficient of the /3-function isn't a complete 
surprise. For an ideal gas, or any system in which the trace of the stress energy tensor 
vanishes, dp/de = 1/3. As is well known, for massless theories the conformal anomaly 
implies that the trace of the stress energy tensor is proportional to the /3-f unction, 
which is why it enters into the relationship between the two viscosities. 

For practical purposes the bulk viscosity is usually neglected relative to the stress 
viscosity. Certainly in the strict perturbative regime ( is negligible relative to rj. 
But the peculiar ratio of integers in front of rj is about 1/5; thus ( is an extemely 
sensitive function of g. As the coupling varies from g ~ 1 to g ~ 3, £ varies from 
essentially zero to a value commensurate with 77. Perhaps this enormous sensitivity 
has phenomenological consequences; after all, g ~ 3 is still a = g 2 /4:ir of order one. 

In summary, while I cannot point to a direct experimental probe of the damping 
rates, they are crucial in sorting out the transport coefficients, and thereby deserving 
of our attention. 

I thank Randy Kobes and Gabor Kunstatter for a most stimulating and enjoyable 
workshop; especially those wonderful lunches on the roof of the Museum of Art! I also 
thank Ben Svetitsky, for correspondence which produced the discussion in sec. I; Rolf 
Baier, for discussions at the workshop which inspired sec. II; and lastly, Eric Braaten, 
with whom the work in sec. Ill was done in collaboration. This work was supported 
in part by the U.S. Department of Energy under contract DE-AC02-76CH00016. 
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